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C ' Abstract 

The cosmological fluid equations are used to study the nonlinear mode coupling of density 
fluctuations. We find that for realistic cosmological spectra there is a significant contribution to 
the nonlinear evolution on scales of interest to large-scale structure from the long-wave part of the 
initial spectrum. A consequence of this mode coupling is that at high redshift, z, the nonlinear scale 
[defined by cr(z) = 1] can be significantly larger than a linear extrapolation would indicate. For the 
OS ' standard CDM spectrum with a as = 1 normalization the mass corresponding to the nonlinear scale 

at z = 20, 10,5 is about 100, 10,3 times (respectively) larger than the linearly extrapolated value. 
We also investigate the possibility of divergent contributions to the density field from long-wave 
modes if the spectral index of the power spectrum n < — 1. Using an approximate non-perturbative 
approach we find that for n > —3 the divergent contribution appears only in the phase. This can 
be related to the large-scale bulk velocity, and clarifies previous results from N-body simulations. 
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1 Power Spectrum at Second Order 



Perturbative studies of the Newtonian cosmological fluid equations have been used to study weakly 
nonlinear effects in the evolution of density fluctuations (e.g., [|J). We have used second-order pertur- 
bation theory to study the enhancement of power and changes in the time-dependence of characteristic 
length scales. We assume that vorticity and pressure are negligible. For the numerical results we have 
used the standard CDM spectrum with Q = 1, Hq = 50kms _1 Mpc -1 and o~g = 1. 
The Fourier transform of the cosmological fluid equtions can be written as: 



+ # = _ d 3 kl t0( kuT )S(k-k U T) , (1) 
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where k is the comoving wavevector, r is conformal time and o(r) [= 1/(1 + z)] is the expansion factor. 
5 and 9 are, respectively, spatial Fourier transforms of 5 = Sp/p and 9 = V -v. Following the standard 
perturbative approach, we expand 5 and 9 in a perturbation series as: 5{k,r) = a8\[k) +a 2 52(k) + .... 
Given the initial power spectrum, and assuming the initial fields to be Gaussian random fields, the 
formal solutions for S(k,r) can be used to obtain its power spectrum to second order: 



P(k,r) =a 2 (r)P 1 (k)+a 4 (T)P 2 (k), 



(3) 



where P\(k) = {8\8\) is the initial (linear) spectrum and P2{k) = (8282 + 28183) is the second-order 
contribution. Pi{k) is an integral over wavevectors involving products of P\(ki) evaluated at different 
k\. We have calculated P(k, r) at different values of a for the standard CDM spectrum and estimated 
the relative contributions from different parts of the initial spectrum. 

We find that the second-order contribution provides a significant enhancement of power in the 
weakly nonlinear regime. It is in excellent agreement with results from a high resolution N-body 
simulation for z J>2, but shows a larger enhancement for smaller z. Due to the effect of long-wave 
mode coupling, the nonlinear evolution responds to the shallower slope of the spectrum at long waves. 
This slows the growth in time of characteristic nonlinear scales, i.e., causes them to be larger at early 
times. We find that the mass corresponding to the scale at which the r.m.s. fluctuation <j(z) = I 
at z = 20, 10, 5 is about 100, 10, 3 times (respectively) larger than indicated by the standard practice 
of extrapolating the nonlinear scale today using the linear spectrum. This result has significant 
implications for nonlinear structures at high z. 

2 Coupling of Long- Wave Modes 

An alternative to the perturbative approach is to approximate the nonlinear terms on the right-hand 
side of the fluid equations shown above so as to include the dominant contribution from the long-wave 
modes. In the limit that the wavevector k\ — ► one finds that the mean square of the nonlinear terms 
takes the form P\{fc)$dk\ P\(k\). This integral is divergent as k\ — > if Pi(fei) oc A;" with n < — 1, even 
while the density variance is finite for n > —3. Indeed some studies of N-body simulations indicate 
that the evolution of n < — 1 spectra is qualitatively different from that of n > —1 spectra, and is 
strongly influenced by the size of the box (the small-/c cutoff) 0, 

The nonlinear terms described above can be approximated by making a Taylor series expansion in 
(ki/k) so as to include the dominant contribution from small k\. On further approximating the small- 
k\ density and velocity fields by the linear fields, the equations reduce to linear partial differential 
equations. We have separated the equations for the phase and amplitude of the density field retaining 
the first order Taylor series term, and have obtained a closed form solution for the phase. 

The divergent nonlinear terms drop out of the equation for the density amplitude implying that 
the power spectrum does not diverge for —3 < n < —1 in this approximation. The solution for the 
phase is divergent for n < — 1. This divergence is related to the divergence of the linear mean square 
bulk velocity, and corresponds to translational motion of the fluid past the fixed coordinate frame. It 
agrees exactly with the deviation of the phases from their initial values found in N-body simulations. 
However this deviation does not seem to indicate dynamically interesting nonlinear evolution as it 
does not affect the amplitude. This conclusion appears to differ from the results of some n = —2 
simulations. We are currently investigating the source of disagreement and also analyzing the long- 
wave limit of perturbative contributions to the power spectrum to verify the above results. 
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